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Abstract 

We consider the random Erdos-Renyi network with enhanced clusteri- 
zation and Ising spins s — ±1 at the network nodes. Mutually linked spins 
interact with energy J. Magnetic properties of the system as dependent 
on the clustering coefficient C are investigated with the Monte Carlo heat 
bath algorithm. For J > the Curie temperature T c increases from 3.9 
to 5.5 when C increases from almost zero to 0.18. Fhese results deviate 
only slightly from the mean field theory. For J < the spin-glass phase 
appears below Tsg\ this temperature decreases with C, on the contrary to 
the mean field calculations. The results are interpreted in terms of social 
systems. 

PACS numbers: 89.65.-s, 64.90.+b 

Keywords: random networks; phase transitions; 

1 Introduction 

Research on systems with random topology was initialized by Paul Flory in 1941 
[UH], but they reappeared only about nine years ago [3] as a subject of great 
interdisciplinary interest with applications in many sciences, from transport to 
biology [U El [E] ■ Most of effort is concentrated on the topological structure 
of growing networks, where the scale-free degree distribution has been found. 
However, progress is achieved in the science of networks as a whole. 

In particular, networks provide a convenient basis of modelling social sys- 
tems [7|. However, the structure of social networks remains debatable. As it 
was discussed by Michael Schnegg [9] , the actual topology of a social network 
depends on social relations in the society. Schnegg investigated six small-scale 
societies, mostly African ones. The degree distribution of three of them were 
found to be close to decreasing functions with some fluctuations, whereas the 
other three displayed a maximum. As the determinant what allows to differ 
between these two groups, Schnegg proposed the amount of reciprocity; when 
it is large, the model degree distribution is close to the Poisson function, known 
to appear in the classical Erdos-Renyi random networks ([HE]. When the reci- 
procity is small, the system is supposed to be close to a scale- free network [5]. 
It seems that the topology of a social structure varies from one society to an- 
other, and various kinds of networks can be appropriate in different cases. We 
note only that, as remarked in Ref. [TU], the clustering coefficient in the social 
networks is much larger, than in the Erdos-Renyi networks. 
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On the other hand, the topology itself does not reflect the richness of be- 
haviour of the social systems, and it is worthwhile to develop a theoretical de- 
scription of these systems with more degrees of freedom. The simplest method 
is to decorate the network nodes i = 1,...,N with additional variables, as 
Ising spins Sj = ±1. These variables are not statistically independent. The 
system total energy is given as E = — J/2 SjSj, where the summation goes 
via all connected pairs of nodes/spins and J is the so-called exchange integral. 
Two basic cases: ferromagnetic interaction (J > 0) which prefers the same sign 
of neighboring spins, and antifcrromagnctic interaction (J < 0), which prefers 
their opposite signs, are to be discussed separately; for a recent review see [IT] . 
In each case, the interaction competed with some noise, commonly modeled as 
the thermal noise. In numerous cases, a variation of the noise intensity leads to 
a phase transition, from the phase with ordered spins (low noise) to the phase 
of disordered spins (large noise) . The amount of noise where the ordered phase 
disappears is comparable to the Curie temperature in ferromagnets. In fact, 
spin degrees of freedom ruled by some stochastic evolution has been used many 
times to describe social processes, as for example the opinion dynamics |12| . In 
these works, ordering of spins can be treated as a demonstration of a collective 
behaviour despite of the presence of noise. The noise level is measured by pa- 
rameter T and usually called temperature. Then the probability of move in the 
system states space from state '1' of energy E\ to the state '2' of energy E2 may 
be given by the probability Pl ^ 2 = [1 + cxp(2A J B/T))]- 1 , where AE = E 2 -E 1 . 
Such dynamic rule is termed as the heat -bath algorithm [13] ■ When we deal 
with the social systems, a society can be modeled by the network and the Curie 
temperature (T c ) can be treated as a measure of an ability of a society to a 
collective action. This ability is expected to depend on the topology of inter- 
personal bonds. In terms of magnetism, the Curie temperature depends on the 
topology of the investigated network [T3] . 

The case of antifcrromagnetic interaction is of special interest for socio- 
physics. As it was argued only recently |15j . a dichotomous behaviour is found 
in some circumstances, as buying or selling |16| . Another example is the dove 
or hawk strategy [TT] . In this case the topology of the social network is par- 
ticularly important because of the geometrical frustration effect, which removes 
the unique and deep minimum of the magnetic energy. Instead, numerous local 
stable or metastable states appear, and the structure of the set of these states 
remains unsolved [18] . We note in passing that the dynamics of the magnetic 
behaviour of antiferromagnetic growing networks was found to be surprisingly 
rich [19]. 

Here we arc interested in the Erdos Renyi random network decorated with 
Ising spins, with ferromagnetic or antifcrromagnctic interaction. To enhance the 
clusterization coefficient, we apply the algorithm proposed by Holme and Kim 
in Rcf. [20] for the growing networks; it is easy to generalize it for the random 
networks. Other methods of network design have been proposed for example in 
Ref. |21j . Our aim is to investigate, how the clustering influences the system 
behaviour. 

The Curie temperature for the random network was obtained already by 
Sumour et al. [22] . However, up to our knowledge the clusterization coefficient 
has not been varied in these networks. Recently, similar problem for J > 
was also considered [221 m the hierarchical network, designed by Watts et al. 
[2~i] . In that network, the clustering was controlled by the probability of linking 
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different clusters. Varying it, one moved from the random network to a set of 
separate clusters. The obtained Curie temperature varied from 3.8J to zero in 
a step- like way, very sharply. Owing to algorithmic/numerical peculiarities of 
the hierarchical network, the size of the system investigated in [53] was limited. 
Also, the clustering led the system to split into subnetworks, what does not occur 
in our case. In the antifcrromagnctic case J < 0, the uncorrelated network is 
believed to show the spin-glass phase at low temperature [57[ HH] • As all J's are 
negative, there is no bond disorder; the frustration is introduced by the disorder 
of the network structure and it is purely geometrical [29] . 



At the first stage of the simulation, a random Erdos-Renyi network is con- 
structed. The starting point is a set of N = 10 6 nodes. A link is placed between 
each two nodes with the probability p. Next, the clustering is enhanced accord- 
ing to the receipt in Rcf. [50] : nodes are selected with two or more neighbours, 
and a new link is placed between each pair of these neighbours with probability 
r selected as to get the average degree (k) = 4. 

The clusterization coefficient C is defined as the average over nodes i = 
1 , . . . , N of the local coefficient d , where 



ki is the degree of i-th node i.e. the number of nodes linked to i, and y, is the 
actual number of links between these ki nodes. The maximal value of Cj is one. 
In our system, the coefficient C varies from almost zero (when p = A/N and 
r = 0) to 0.18. This enhancement appears to influence the degree distribution. 
As it is shown in Fig. Q] the plot is Poissonian for small C, but deviates from 
this curve for C = 0.18. 

The magnetic behaviour of the system is also influenced by the clusterization. 
Let us start from the ferromagnetic case, where the exchange integral J is set to 
+1. We calculate the average magnetization M = Si/N against temperature, 
M(T), for r = 0. Initially, all spins are set to +1, and the time average is taken 
after N t time steps from another N t time steps. One timestep is equivalent 
to N attempts to flip a spin. For the comparison with the results of [55], the 
calculations for (k) = 4 are repeated for (k) = 3 and (k) = 2, N = 10 4 . The 
obtained curves for the squared magnetization are presented in Fig. [5] As we 
see, the obtained Curie temperature increases with (k). Also, the mean-field- 
character of the system is confirmed, as it was concluded in [22], because the 
curves M 2 (T) are approximately linear. For (k) = 2, the value of the Curie 
temperature, T c /J « 1.7 agrees with the one obtained in [55J. For (k) = 4, 
the value of the Curie temperature T c /J m 3.8 agrees with the result of [53] 
for the case of random network, where the homophily parameter a = — In 2 
(for the discussion of the interpretation of a see Ref. [24]). We note also 
that the analytical solution [T3] for the uncorrelated random networks with the 
Poissonian degree distribution 
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Figure 1: The degree distribution for the random network with C ~ (stars) and 
C = 0.18 (circles), for TV = 10 6 nodes. The former is the Poisson distribution 
for (k) = 4 and C « 0; there, <fc 2 } = 20.0 and z 2 = 16.0. For C = 0.18 the 
simulation gives (k 2 ) = 25.5 and z 2 = 21.5. 



gives T c /J = 3.91, 2.88 and 1.82 for (k) = 4, 3 and 2, respectively. The obtained 
values are slightly smaller. The differences are due mainly to the finite-size 
effect. 

These calculations were repeated again for (k) = 4, J > 0, r > and a 
larger network, i.e. TV = 10 5 . The results shown in Fig. [3] indicate, that the 
Curie temperature increases with the clustering coefficient C . This increase does 
follow from the increase of the second moment (k 2 ) of the degree distribution. 
As shown in Fig. [H the theoretical estimation of Eq. differs only slightly 
from the observed value of T c . The same is true for the more general formula 

m 

i = \^(P^f) (3) 



T c 2 \z 2 - zi 

where Z\ = (k) = 4, and z 2 is the number of second neighbours. Both equations 
2 and 3 give the same theoretical values. Indeed, our numerical calculations 
confirm the relation z 2 ~< k(k — 1) >. 

In the antifcrromagnctic case (J < 0) the disordered topology of the random 
network does not allow to investigate the staggered magnetization as the order 
parameter. We made an attempt to calculate the magnetic specific heat Cy 
against temperature. For the lattice antifcrromagnets, it is the maximum of Cy 
at given temperature what marks an existence of the phase transition. Indeed, 
we found an indication of such a transition for the case r = 0, where the clus- 
terization coefficient vanishes. Slightly below the maximum, Cy obtained from 
the derivative of energy with respect to T starts to differ from Cy calculated 
from the variance of energy. This split is due to the lack of thermal equilibrium 
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Figure 2: Thermal dependence of the squared magnetization for J > 0, N = 10 4 
nodes, r = 0. The mean node degree (k) is 2 (left curve), 3 (middle) and 4 (right 
curve). (iViter = 2-10 4 ). 



in the spin-glass phase. Still, the position of the maximum can be determined. 
In Fig. [5j only the result from the thermal derivative of energy is shown. These 
calculations are preceded by N t = 10 4 time steps, till the time dependence of 
energy reaches the plateau; this was checked to be long enough if T > 0.7. 

The plots shown in Fig. [5] show that the transition temperature T$g de- 
creases with the clustering coefficient C. This result is confirmed by the calcu- 
lations of the Edwards- Anderson spin-glass order parameter q [25l [26] 

s = ^£(;X>«>) 2 w 

i t=l 

where Si(t) is i-th spin at time t. The obtained plots for q(T) for various C are 
shown in Fig. [51 Here also, the temperature where q vanishes decreases from 
1.9 for C — to about 1.4 for larger values of C. 



3 Discussion 

New results obtained for the ferromagnetic case indicate that the enhancement 
of the clusterization by additional links does not invalidate the assumption of 
the lack of correlations, which is usually made in the theory [14] of random 
networks. However, for C = 0.14 and above this value, the loops introduced to 
the system slightly enhance the Curie temperature, as shown in Fig. 4. 

For the antiferromagnetic case, the situation is more complex. For negligible 
clusterization, the maximum of the specific heat indicates the phase transition 
between the paramagnetic and the spin-glass phase. The transition temperature 
Tsg r °ugly accords with the mean field result T$g = 1-91 from Eq.5. 
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Figure 3: Thermal dependence of the squared magnetization for J > 0, different 
values of C, and N = 10 5 nodes. The mean node degree (k) = 4. The clustering 
coefficient C « 0, 0.05, 0.09, 0.14 and 0.18 from left to right. 
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Figure 4: The Curie temperature obtained from the simulation (crosses) and 
from Eq. © or © (circles). 
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Figure 5: The magnetic specific heat Cy(T) for the case J < 0, where C—0, 
0.09 and 0.18 (crosses, circles and triangles, respectively). 
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Figure 6: The Edwards- Anderson order parameter q(T) for the case J < 0, 
where C=0, 0.09 and 0.18 (crosses, circles and triangles, respectively). 
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which follows from Eq. 92 in Rcf. [TT] and from our delta-like distribution 
p(J) = S(J + 1) of the exchange integral J. However, the theoretical values 
of T SG increase with C to 1.92 for C = 0.09 and to 2.17 for C = 0.18. Our 
numerical results on Cy and q indicate, that actually T$g decreases with C. 

The topology of the system allows to treat it as a disordered version of 
the Archimedean lattice, discussed in Ref. [55]. There, the triangles of spins 
were connected as to form the periodic structure. The ground state degeneracy 
of this Archimedean lattice increases exponentially with the system size. The 
same rule applies to our disordered lattice if the density of the triangles is large 
enough. 

As discussed in Ref. [TT], there are links from the problem of spin-glass in 
disordered networks to some interested NP-complete problems as MAX-CUT or 
the satisfiability. It seems that some counterparts of our system can be found 
also in the game theory, and in particular in the network congestion games 
[30] : for a recent review of games on networks see Ref. [31]. To maintain our 
social interpretation, we can state that the case of a ferromagnet the smooth 
increase of the Curie temperature with the coefficient C marks the fact that the 
ordered phase is less fragile if the system has the clustered structure. Consid- 
ering an antiferromagnet we take into account that effects of social interactions 
are sometimes incoherent, that the society can be polarized in some way and 
that it is possible that increasing the density of social ties can be in conflict 
with an ability to find best solutions to everyday problems. In this sense, geo- 
metrical frustration combined with noise can reproduce some social phenomena 
which disable a coherent social action. A next step to reflect features of the 
society should be to consider a distribution of positive and negative interaction 
of different sizes. This, however, requires much larger conceptual and numerical 
effort. 
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